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SECTION I 


INTRODUCTION 


There are spectral regions in the infrared that are not cccluded by 


atmospheric absorption’. The coincidence of good detectors and atmospheric 
transparency in the 2-6 ym and the 8-12 um bands makes these .pectral regions 
of special interest. Military applications in the 2-6 um band and, to a 
lesser extent, the 8-12 um band have been discussed elsewhere”, Prior tu 
1970, tunable sources had not been developed to provide continuous coverage 
of these bands, although discrete line lasers were available at fixed wave- 


lengths throughout the bands. Sauce then, the availability of new nonlinear 
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optical materials with wide band infrared transparency has mcde possible the 
construction of cptical parametric oscillator (OPO) and mixer devices that 
operate contiruously over the 2-6 ym and &-12 um bands. (Although tunable 
semiconductor sources such as the lead sait diodes and spin flip Raman lasers 
have recently veen developed, we do not consider such devices in this repert.) ‘ 
In Section II, we discuss parametric oscillators. A review of the chrono- 7 
logica) development of wavelength coverage and ¢ discussion of power output 
and efficiency vs wavelength are given. In Section III, we review recent work 


on turable mixers. In Section IV, we calculate and discuss the theoretical 


performance of 11 nonlinear optical material. In Section V we show, as an 
example, how the theoretical performance data can be used in selecting a 
pump wavelength and a nonlinear optics. naterial for a singly resonant OPO 


operating in the 3-6 um band. In the Appendices we include the refractive 
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index data and Sellmeier equations for each nonlinear material and present 


vo med 


the details of the co:sputer program used for the calculations. 
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PARAMETRIC OSCILLATORS | 


A parametric os: illator generates two waves, a signs] and an idler, 


from noise in a nonlinear crystal when properly excited by a leser pump 


j 
beam, The energy and momentum conservation requirements on the pump, | 
signal, and idler beams are: j 

w, sus tw; Energy Conservation (1) 
ke ket k Momentum Conservation (2) 


The theoretical maximum efficiency of the process for generating signal 


energy is w/a, and for idler energy ic w/e, (Figure 1). For the 
degenerate case, t.e. “su w/?, 100% conversion is possible to a 
single line, In practice, 67% conversion has been obtained in a degenerate 
oscillator® and 100% conversion has been reported in the case of degenerate 
upconversion where the frequency of the pump is doubled. * Detailed dis- 
cussions of parametric oscillators are contained in Ref orences 5 and 6. 

The cnronological progress in spectral coverage for nonlinear devices 
is shown in Figure 2 for the visible and infrared region, The numbers in 
Figure 2 correspond to the reterence numbers given in Table I and in the 
list of references at the end of this report. LiNbO, is the most commonly 


exploited material and has been used from the green to 3.7 um where it is 


limited by absorption. The materials used for parametric oscillators must 


; be transparent or low loss at the pump wavelength as well as the signal and 


idler; they must also have adequate optical quality to allow an oscillator 
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cavity to be formed. Pumping LiNbO, with different visible wavelengths 
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Theoretical Conversion Efficiency as a Function of the 


Figure 1, 
Ratio of the Generated Frequency to the Pump Frequency 
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yielded spectral coverage from 0.548 to 3.65 um; more recently, pumping 


LiNbO, with an electronically tunable dye laser™” offered rapid selection 


of wavelengths in the near IR. Proustite (Ag,AsS,) pumped by Nd: Cawo, 
* 
provided coverage out to 8.5 um; 28 CdSe pumped by CaF, :hy** reached 


30. Ammann, Yarborough and Falk?? achieved simultaneous optical 


13.7 um 
parametric oscillation and frequency doubling of the signal in LiNbO, pumped 
by 1.06 um to generate wavelengths higher than the initial pump. Hence, the 
infrared and visible spectrum of interest can be covered by parametric 
oscillators. 

The highest IR peak powers (Figure 3) were obtained from parametric 


osciliators pumped by rupy lasers, 11.15,15,16 


pumped by CaF, :Dy**,5° Peak powers up to 15 kilowatts have been obtained 
from parametric oscillators pumped by neodymium lasers. In contrast, all 
the h:gh average power parametric oscillators have been pumped by neodymium 
lasers capable of high repetition rate performance (Figure 4). The only - 
parametric escillator shown in Figure 4 that was not pumped by iiecdgruitia Was 
a CW argon laser-punped lithium niobate oscillator. 4 It should also be 
noted that the parametric oscillators of References 3 and 20 were internal 
parametric oscillators that exposed the oscillator to the total laser flux 
in the cavity. The parametric oscillator of Reference 19 was pumped in 

the greer. which resulted in 90° phase matching for high average sewer 
performance, 

In Figure S$, power conversion efficiency for the various parametric 
oscillators is shown. The highest efficiency was that or the previously 
mentioned degenerate oscillator” at 2.16 wm. Power conversion efficiency of 
46% from a 0.659 um pump to two tunable outputs was reported. 2? A 4st peak 


power conversion to the signal was obtained in a ruby-pumped singly resonant 


oscillator, !® Conversion efficiency of 40% as measured by pump depletion 
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followed by the CdSe oscillator 
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Figure 3. 
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Peak Power vs Wavelength 
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Figure 5 
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was reported?” 


for an oscillator with the signal resonated. The various means 
of reporting power conversion efficiencies were not totally self consistent, 
and in some cases, an interpretation was required. For explicit details on 
the efficiencies see the references and Table I. 

There are three ways commonly cited to give an idea of the power 
conversion efficiency, Average power conversion efficiency is the ratio of - 
the average output power of the parametric oscillator to the average input 
power of the pump. Peak power conversion refers to the ratio of the peak 
power of the parametric oscillator to that of the pump. Note from Table I 
that the pulse length of the oscillator is frequently shorter than thet of 
the laser pump; this results in an apparent higher efficiency for peak power 
conversion because of pulse compression. Pump depletion is a measure of 
the amount of energy lost from the pump in passing through the parametric 
oscillator. This lost energy is converted to signal and idler energy in 
the parametric oscillator as well as absorbed and scattered in the parametric 
oscillator elements. Hence, not all energy measured by pump depletion is 
transformed to useful parametric oscillator output. 

Parametric oscillators of lithium niobate have been tuned by angle, by 
temperature, and by the electrooptic effect (Table I). The bringing of 


the crystal to thermal equilibrium at new temperatures changes the indices 


of refraction, and hence, the frequency of the parametric oscillator. This 


process is a slower tuning technique than that of merely rotating the ( 


crystal. The technique of changing parametric oscillator frequency by 
changing the wavelength of the dye laser pump has been recently demonstrated 
using an acousto-optic filter, °° This offers random frequency selection 


within the limits of the dye laser pump and maximum flexibility since 


10 


frequencies can be shifted on a pulse to pulse basis. FElectrooptic effect 
tuning! demonstrated a frequency shift of 6.2 A°/kV/cm in LiNbO,, so it 
cannot be effectively used for large frequency excursions because of the 
voltage involved. 


2? “The advantage of a singly resonant 


Bjorkholm succinctly puints out 
oscillator is thet it should be much more nearly continuously tunable than 
a doubly resonant oscillator. In the latter, because of the accidental 
distribution of the signal and idler cavity modes, frequency excursions 
from the operating point predicted by index-matching data can be large 
with respect to the axial mode spacing of che oscillator (~0.1A°). For 
the singly resonant oscillator, the tuning discontinuities should te no 
larger than one-half the mode spacing." In another article!® Bjorkholm 


remarks that the tunability of the SRO is at least five times and probably 


an order of magnitude more accurate than the tunability of the DRO. 


SECTION IIT 


TUNABLE MIXERS 


Coherent infrared radiation can also be produced as the difference fre- 


quency generated when two coherent sources, the pump and the signal waves, 
are mixed in a suitable nonlinear crystal at the proper phases matching 
angle. Table .I lists typical pump wavelengths and correspor.ding signal 


wavelengths required to generate IR wavelengths from 2 to 12 um. Again, 


the energy and momentum relationships (, - us say and x k, = Lo 


respectively) must be satisfied to have efficient power transfer. 


Most experiments to date have utilized single pass systems which do 


ee et ee 


not require vptical cavities and associated mirrors. Such systems offer 
the advantage of simplicity, but sacrifice the more efficient utilization 
of input power through multiple reflections in the cavity. However, crystals 
‘of very high optical quality are required for use in optical] parametric 
oscillators, whereas poorer quality crystals can be used in. wixers. 

In January 1971, Dewey and Hocker! reported infrared difference fre- 
quency generation by mixing the inputs of a Q-switched ruby laser and a 
tunable dye laser in a LiNbO, crystal, In their experiment, the output of 


: an unpolarized Q-switched ruby laser was separated into two parts by a beam 


Splitter. One beam pumped a tunable dye laser. The output from the dye 


laser was then combined with the second part of the ruby beam and directed 


into the LiNbO, crystal which was rotated to maintain the proper phase match- 


: 
4 
y 
| 
! 
1 


ing angle for different wavelengths. Pulses of approximately 6 KW of peak 


power were produced in the ~ to 4 um range for an initial input of 4 Mw from 


the ruby laser. 
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The dye laser mixing technique waz extended further into the If regiu.: 


by Hanna, Saith, ano Stanley’~ whon they reported (December 1971) that they 
had mixed the outputs of a ruby Jaser and a tunable dye laser (pumped by 
the ruby laser) in a proustite crysta! to generate peak powers of a few 
watts at approximately 5 um and over.a range of 10.1 to 12.7 um. 


In June 1972, Meltzer and Goldberg”? 


obtained tunable difference fre- 
quency generation in Lil0, over a range of 4,1 to S.2 um. They mounted the 
Lilo, crystal inside the dye laser cavity. The dye cell was pumped by a 

ruby leser and tuned with a diffraction grating. Phase matching was achieved 
by rotation of the crystal. IR output pulses of spproximately 100 W peak 
power at 4.7 um were obtained for a total ruby laser input of 4 MW. The dye 
laser power was approximately 70 kW. 


i reported the 


In December 1972, Bhar, Hanna, Luther-Davies, and Smith 
use of an GPC as their tunable source to generate two input wavelengths 
which, when mixed in a proustite crystal, produced tunable IR radiation over 
a range of 7.8 to 11.9 um. The OPO was tuned from 1.87 to 2.47 ym. The 
total input power to the crystal was approximately 100 W. 

Down-converted IR radiation from 4.6 to 12 um was reported by Hanna, 
Rampal and Smith in June 197% when they mixed a ruby laser beam in silver 


thiogallate (AgGaS,) with the output of various ruby-pumped dye lasers. The 


Q-switched ruby laser produced a peak output power of 1 MW (TEM, mode) with 
a pulse duration of 10 nsec (fwhm). Part of this power was separated by a 
beam splitter and was used to pump the dye laser. Four separate dyes were \ 


used to cover the 4.6 to 12 um spectral range. Peak powers of hundreds of ! 


milliwatts were produced. 


\ 
4 
14 | 
| 


Decker : nd Titter??9°5 reported (Apris 1973 and July 1973) tunable 
difference frequency generation in proustite over the range 3.2 to 5.65 


wa with a Q-switched rub, lases and tunable dye jaser and over the range 


oeeudee 


§,.82 to 7.25 wa with the outputs of two independent ruby pumped dye lasors. 
Peak IR powers up to several KW were obtained. 

Byer, Choy, Herbst, Chemla, and Feigelson”” recently reported (January 
1974) continuously tunable IR difrerence frequency generation in silver 


gallium selenide (AgGaSe,) over the range 7 to 13 wm by mixing the outputs 


from a LiNbO, parametric osciliator. A 1.318 wm output from a Nd:YAG laser 


was doubled by a Lil0, crystal to obtain 0.659 um which in turn puaped a 


ela salto 


LiNbO, temperature-tuned oscillator to produce a 1,62 to 1.47 um signal. 


- ROD Me Relate ele enr me aE nee mn an Hy ee ce eg 


This signal was matched with the 1.318 um pump to produce the 7 to 13 um 


9 


IR radiation. 

Tunable coherent infrared radiation produced in a four-wave mixiig 
process was reported in April 1973 by Sorokin, Wynne, and Lankard’’, The 
mixing process was achieved by utilizing resonant enhancement of a third 
order optical nonlinearity in an alkali metal vapor (potassium at 500°C). 
A 100 «W nitrogen laser was used to simultaneously pump two dye lasers 


whose outputs were focused into a heat pipe oven containing the metal 


vapor, The beam from the first dye laser at frequency w,, which was tuned 


a. 2 saa Cath NN, cr at iD ees. 


to the 4s-Sp resonance lines of potassium (0.404 um), induced radiation (a 


Stokes wave) at frequency W> by stimulated Raman scattering between the 


ee a gett 


Sp-Ss levels. The beam from the second dye laser at frequency uw, then 


0 ae 


“beat" together with the first two waves to create a fourth wave at fre- 
quency w, (in the infrared) such that Wy Ty = Wy - wy (the energy matching 


relationship). The corresponding momentum relationship that had to be 


Led putida tee 


1S 


Sede 


satisfied was dnjw, © 4njw, -Anow-A2,05 where Ar. is the dispersion in the 
refractive index. By changing the lesex dyes and by mixing small amour.ts 
of other gases with the potassium vapor, it was possible to obtain continu- 
ously variable wavelengths from 2 to 20 um. IR peak power outputs up to 100 
WN were obtained. 

Mixing experiments reported to date have demonstrated the feasibility 
of generating tunable infrared coherent radiation at very low power levels 
in a number of differe;t materials. For the most part, pump sources have 
been single. ahot or low repetition rate lasers with outputs in the visible. 
To achieve worthwhile power levels requires that the two input sources, one 
of which must be tunable, operate in the 1-2.5 um region for increased 
efficiency. The lack of a suitable high power tunable source at these 


wavelengths is a serious limitation at present. 
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SECTION IV 
THEORETICAL PERFORMANCE OF NONLINEAR OPTICAL MATERIALS 


The theory of parametric interaction of confocally focused light beans>* 


was used to provide a basis for comparing materials for parametric generation 
of tunable radiation in the 2-6 and 8-12 um regions. The effects of walkoff 
due to birefringence were taken into account. 

The materials used in the comparison are listed in Table III along with 
the nonlinear coefficients used in the calculations. The refractive indices 
used were obtained from the references also listed in Table III. For com- 
pleteness we have coapiled in Appendix A the room temperature index data and 
Sellmeier equations for each material. The list of materials is not com- 
prehensive but does contain materials already used in nonlinear optical devices 
as well as several promising materials recently synthesized. Of the 11 
materials listed in Table III, only LiNbO, and LiI0, are commercially available. 
All the other materials are in various stages of development. Of the nine 


non-commercially available materials listed in Table III, Loth Ag,Ass,'6+24,28 


ana CdSe20» 27.50 


have been used in infrared parametric oscillators. Also, 
33,3 3 
sA8S5,°1°59995 agcas., 


These materials, then, may be considered to be the most advanced of the nine 


Ag,AsS . and AgGase,”” have been used for infrared mixing. 
experimental nonlinear materials, 

In Figures 6 to 12, curves of the parametric gain per watt of pump power 
as a function of oscillator output wavelength are given for the pump wave- 
lengths 0.694, 0.946, 1.06, 1.3, 1.83, 2.1, and 5.3 um, respectively. All 
of the materials of Table III that can be phase matched and which are trans- 


parent at a given pump wavelength are included in each figure. A crystal 


length of one cm was assumed in all cases. Each curve extends over 
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TABLE III NONLINEAR OPTICAL MATERIALS 


Material] Crystal] Birefrin- 
Class 
YP 


NLO Coefficients 
(x 10*/e, m/V) 


Refractive 
Indices 


ld, = 


lay! = lage 


12 (Ref 39) | Ref. 39,40 


AgGaSe., 42m | Negative Id, = 33 (Ref 41) | Ref. 41,42 
CdGeAs. 42m | Positive | [d,,] = 236(Ref 43) 
HgS 32 Positive Id,,I = 5U (Ref 46) | Ref. 47,48 
Ag.ASS., Negative 25.6(Ref 49) | Ref. 48,50 
lana erie 
Ag ,SbS,, Negative = 13.4(Ref 51)} Ref. 52 
= 12.6 
LiNbO, Negative 3.6 (Rbf 53) | Ref. 48,54 
id ing oe 
T1,AsSe, Negative = 44.2(Ref 55) | Ref. 55 
Sa inna 0 
Lil0, i |a,,| = 7-5 (Ref 56)| Ref. 56,57 
. Ic), = 0 
CdSe 6mm "| Negative Id, = 19 (Ref 58) | Ref. 31 
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Figure 6, Parametric Gain Per Watt vs Output Wavelength for a Pump 
Wavelength of 0.694 um. A Crystal Thickness of 1 cm is 
Assumed. Curves are for Type I Phase Matching Unless 
Otherwise Noted. 
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Figure 8. Parametric Gain Per Watt vs Cutput Wavelength for a Pump 
Wavelength of 1,06 um. Curves are for Type I Phase 
Match.ng Unless Otherwise Noted. 
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Figure 9. Parametric Gain Per Watt vs Output Wavelength for a Pump 
Wavelength of 1.3 um. Curves are for Type I Phase 
Matching lnless Otherwise Noted. 
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Figure 10, Parametric Gain Per Watt vs Output Wavelength for a Pwap 
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Figure 11, Parametric Gain Per Watt vs Output Wavelength for a Pump 
Wavelength of 2.1 um. Curves are for Type I Phase i 
Matching Unless Otherwise Noted. 
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Figure 12, Parametric Gain Per Watt vs Output Wavelength for a Pump 
Wavelength of 5.3 ua. Curves are for Type I Phase 
Matching Unless Otherwise Noted. 
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the transparency range of the material. in the case of the phosphides and 
sulfides, however, the longest prsctical wavelength is limited by multiphonon 
absorption to about 10 wa. The Roman numerals beside each msterial identify 
the type of phase astching (see the dofinitions given in Table IV). In 

cases where both Type I and Type II phase matching are possible, only the 
higher gain curve is shown if (as, for example, in the case of Ag,AsS.) the 
gains for the two types differ by a large factor. In cases where the gains 
for both types of phase matching are comparable, both curves are given. In 
previously published comparisons of materials,>’° parametric gains are dis- 
Played as horizontal lines vs wavelength (deroting only the transparency 
renge), and different pump wavelengths are included on the same diagram. The 
actual gain curves, their crossings, and their relative shifts froa one pump 
wavelength to another, as seen in Figures 6 to 12, show that a simple composite 
diagram can be misleading. 

High gain implies a low pumping power required to reach oscillator 

threshold. Thus, to build a practical device, one must select a material 

with sufficient gain to reach threshold and, ideally, saturation at a power 
level within the state-of-the-art of the pump lasor. However, gain is a 
necessary but not a sufficient criterion on which to base the choice. If the 
high gain can be achieved only dy such tight focusing that the pump power 
density exceeds the material's damage threshold, then obviously the material 
cannot be used, Other criteria are availability and optical quality of the 
material, and the availability of suitable anti-reflecting and high reflecting 
coatings. | 

To answar the question of power density requirements, we present a de- 

tailed comparison of materials at the pump wavelengths 0,694, 1.06, 1.3, and 
2.1 um in Tables V to VIII, respectively. These pump wavelengths are represent - 
ative of the ruby, na", and Ho* solid state lasers. 
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TABLE IV DBFINITION OF PHASE MATCHING TYPES 


Beam Polarizations 


Birefringence Phase Matching 
Type Type 
Signal Idler Pump 
Q, > dp/2) 


Positive (n° > n°) e ° 
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COMPARISON OF PARAMETRIC OSCILLATOR PERFORMANCE OF MATERIALS AT A PUMP WAVELENGTE OF 9.694 2 
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COMPARISON OF PARAMETRIC OSCILLATOR PERFORMANCE OF MATERIALS AT A PUMP WAVELENGTH OF 1.06 us 
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VII. 
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OF MATERIALS AT A PUMP WAVELENGTH OF 1.3yum 
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The first coluan in Tables V to VIII identifies the material and the 
type of phase matching. The second and third columns give the output wave- 
lengths, The fourth column i. che gain per watt for a 1 cm thick crystal. 
The next two columns give the pump power required to reach threshold for 
singly and doubly resonant oscillators according to the following relation- 


ships:> 


(gain/watt) x Pp" 2a (SRO) (3) 


s 
(gain/watt) xX Pe = 9,0; (DRO) (4) 


where a. and a, are the single pass power losses at the signal and idler 
wavelengths, respectively. For the comparison presented in Tables V to VIII, 
we have assumed values of asa, 5 0.1. At degeneracy Q, = Av), a singly 
resonant oscillator has the same threshold as a doubly resonant oscillator. 
Hence, the numbers given in parentheses in the SRO columns are fictitious 
and should be interpreted as the pump power for an SRO operating near 
degeneracy. 


The next two columns give the pump power required for pump depletion 


(saturation) according to” 
2 
PS = (1/2) P, (SRO) (S) 
“en = 4P. (DRO) (6) 


In Equations S and 6, it is assumed that the pump pulse length is much 
greater than the time required to build up the signal from the incoherent 
background noise, In cases where buildup time is not negligible, the pump 
powers required for high conversion efficiency depend on the gain per ass 
and the ratio of pulse length to oscillator cavity round trip time .° Since 


these parameters can vary over wide ranges, the problem of buildup time is 


ta Pate |e 
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not treated here. However, a good rule of thumb is that a parametric 
oscillator pumped by pulses shorter than -100 ns may be expected to be 
limited by buildup time. 

The columns labeled A in Tables V to VIII cre the cross-section areas 


of the Gaussian beams at the cavity confocal point and are given by 
2 
Ae wi /2 (7) 
° 
= bor! (nyen,) 


where, for a 1 cm crystal, 


b= B7/2. at > 2 (8) 
2 
do, al B < 2 
p? mon, /4h, (9) 


In these equations, Wy is the Gaussian beam waist, » is the pump wavelength, 
ny» 2, and ny, are the refractive indices at the signal, idler, and pump 
wavelengths, and p is the walkoff angle. Equatians 8 maximize the confocai 
parameter b (and, hence, minimize the power density at the focus) for a 

given value of the walkoff parameter B. Values of bo larger than given by 
Equations 8 are, to a good approximation, offset by a reduction in gain, so the 
cross-sections in the Tables V to VIII are essentially optimum. 


The columns following A give the power densities at threshold and at 


saturation. The final two columns give the theoretical maximum conversion 


efficiency at the two output wavelengths. 
f 
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SECTION V 
RANKING OF MATERIALS POR A SPECIFIC APPLICATION 


As an example of how the detailed data presented in Section IV can be 
used, let us consider a specific application, namely, a parametric oscillator 
with output in the 3-6 um region, at a prf of S kHz in an SRO configuration. 


To achieve this repetition rate at the power levels needed to drive a para- 


metric oscillator requires either a na>* or Ho>* pump laser operating in a 


Q-switched, cw-pumped mode. To obtain the highest possible average power, 


_ laser should operate on the 1.06 um (1.08 wm) line of YAG:Nd>* 


3¢ 


the Nd 
(YA10,:Nd5*), rather than at other possible Nd wavelengths (0.94 ym, 1.3 um, 
etc.). At an average power output of 10 watts TEM, (a reasonable upper limit 
for the present state-of-the-art), a train of 100 nsec pulses at 5 kHz will 

have a peak power per pulse on the order of 20 kW at the 1.06 um or 1.08 um 
pump wavelength. The lio” * laser operates at 2.1 wn. A reasonable upper limit 
for state-of-the-- 7% 1AGia-He laser output is 2 watts TEM 0° Thus, the pump 
limit is about 4 kW per pulse at the 2.1] um pum wavelength. 

Using the data in Tables VI and VIII, we can rank the parametric oscillator 
materials in the order of their merit for the application just posed. This 
ranking is displayed in Tables IX and X. For purposes of ranking materials in 
Tables IX and X, we have assu: dam, ‘hreshold power densities of 100 MW/cm? 
for LiNbO, and Lilo, and 50 MW/cm? for all other materials. 

According to Tables IX and X, HgS and ZnGeP., are the best materials for 
pumping at 1.06 and 2.1 um, respective) Both materials have peak power and 
power density requirements well below the upper limits discussed previously. 


For SRO configurations, AgGaSe, pumped at 2.1 ym is the only other material 
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TABLE IX 


LIMITATIONS OF MATERIALS FOR 1.06 ym-PUMPED SINGLY RESONANT PARAMETRIC OSCILLATORS 
OPERATING IN THo 2-6 ym OUTPUT RANGE AT S KHz PRF 


Material 


Hgs 


Threshold Power 
Less than 20KW? 
Yes 
Yes 
Yes 
Yes 
Yes 


No 


* but not beyond © Spa 


b 


but not beyond ~ 4ym 


Saturation Power 
Less Than 20 KW? 


Threshold Power 
Density Less Than 
Damage Threshold? 

Yes 

Yes 

Yes 

Yes 

Yes 


No 


© limited by absorption to wavelengths shorter than » 4.5m 
q but not beyond © Sum 
© limited by absorption to wavelengths shorter than +. 5.5m 


ete nea een tne en 
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Saturation Power 
Density Less Than 
Damage Threshold? 


Yes 
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TABLE X 


LIMITATIONS OF MATBRIALS FOR 2.1 ys-PUMPBD SINGLY RESONANT PARAMETRIC OSCILLATORS 
OPERATING IN THE 3-6 ym OUTPUT RANGE AT S KHz PRP 


Threshold Power Saturation Power 
Threshold Power Saturation Power Density Less Than Density Less Than 
Material Less Than 4 KW? Less Than 4 KW? Damage Threshold? Damage Threshold? 


zinGeP. Yes Yes Yes Yes 
AgGaSe. Yes Yes Yey Yes 
HgS Yes No Yes Yes 
T1,AsSe, Yes No Yes Yes 
Ag,AsS, No No Yes Yes 
Ag ,SbS,, No No Yes No 
AgGaS, No No No No 
I 
t 
k ; 
‘ 
41 
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that has saturation power and power density levels within the state-of-the-art. 
However, most other materials can.be used in SRO configurations up to some 
level between threshold and complete pump depletion. The materials that can 
be categorically excluded from consideration in SRO configurations are Lilo, 


pumped at 1.06 um and AgGaS,,, Ag,SbS,, and Ag,AsS, pumped at 2.1 un. 
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SECTION VI 
SUMMARY AND CONCLUSIONS 


In this report, we have reviewed the state-of-the-art of parametric 
oscillators and mixers, and have provided detailed theoretical performance 
gata on the most promising materials for use in the 2 to 6 and 8 to 12 um 
regions, and have performed a critical evaluation of the materials usable 
in the 2 to 6 um region. The spectral coverage available from parametric 
oscillators and mixers has increased dramatically since 1970. This has 


been made possible by the development of new nonlinear optical materials 


lhl et cae at, sae nal ln ain alt len athe 


such as Ag,ASS,, CdSe, and AgGaSe.. High average powers and conversion 
efficiencies have yet to be obtained over the entire infrared bands of 
interest. Substantial improvements in power and efficiency can be made by 


further development and use of materials such as HgS and ZnGeP,. 


43 


Wien itary: Galata EERE ce 7 ERB A wh Sin Sone tis ERAT Neat TOONS i Ne} ts, teat IAN rae eect Rant cea la nce eee me 


10, 


REFERENCES 


W.C. Eppers, Jr., "Atmospheric Transmission,"' Ch. 3 of Robert J. Pressley, 
CRC Handbook of Lasers with Selected Data on Optical Technology, 

Chemical Rubber Company, (1971). 

Proceedings, Sixth DOD Conference on Laser Technology, March 1974 (to be 
published). 

B.0. Ammann, Three to S Micron Parametric Oscillator, AFAL-TR-72-210, 

May 1972. 

J.B. Geusic, H.J. Levinstein, S. Singh, ?.G. Smith, and L.G. Van Uitert, 
“Continuous 0.532-ym Solid-State Source Using Ba2NaNb-0,¢, " Appl, Phys. 
Letters 12, 306 (1968). 

S.E. Harris, ‘Tunable Optical Parametric Oscillators," Proc. IEEE 57, 
2096 (1969). 

R.L.e Byer, "Optical Parametric Oscillators," Stanford Microwave Laboratory 
Report 2140, January 1973. Published in Treatise in Quantum Electronics 
(Herbert Rabin and C,L. Tang, eds), Academic Press 1973. 

J.A. Giordmaine and R.C. Miller, "Visible Coherent Parametric Oscillation 
in LiNbO, at Optical Frequencies," Phys. Rev. Letters 14, 973 (1965). 
S.A. Akhmanov, A.I. Kovrigin, V.A. Kolosov, A.S. Piskarskas, V.V. Fadeev, 
and R.V. Khokhlov, "Tunable Parametric Light Generator with KDP Crystal," 
J. Exptl, Theoret. Phys. Letters 3, 241 (1966). 

J.A. Giordmaine anc &.C. Miller, "Optical Parametric Oscillation in the 


Visible Spectrum," Appl. Phys. Letters 9, 298 (1966). 


R.C. Miller and W. A. Nordland, ''Tunable LiNbO, Optical Parametric Oscillator 


with External Mirrors,'' Appl. Phys. Letters 10, 53 (1967). 


44 


! 
i 


vine 


Mn 


CES nee 


21. 


22. 


a alti El a te tte Ml a ta eee teh eee 2 


A ma. 
ee a TT a a ait 


L.B. Kreuzer, "Ruby-Laser-Pumped Optical Parametric Oscillator with 
Blectro-Optic Effect Tuning," Appl. Phys, Letters 10, 336 (1967). 
R.G, Smith, J.6. Geusic, H.J. Levinstein, J.J. Rubin, S. Singh, and 
L.G. Van Uitert, "Continucus Optical Parametric Oscillation in 
Ba_NaNb,0,,," Appl. Phys. Letters 12, 308 (1968). 
J.E. Bjorkholm, "Bfficient: Parametric Oscillation Using Doubly and Singly 
Resonant Cevities," Appl. Phys. Letters 13, 53 (1968). 
R.L. Byer, M.K. Oshman, J.F. Young, and S.E. Harris, "Visible CW Parametric 
Oscillator," Appl. Phys. Letters 13, 109 (1968). 
J.E. Bjorkholm, "Spectral Properties of Doubly and Singly Resonant Pulsed 
Optical Parametric Oscillators," Appl. Phys. Letters 13, 399 (1968). 
J, Falk and J.E. Murray, "Single Cavity Noncollinear Optical Parametric 
Oscillation," Appl. Phys. Letters 14, 245 (1969). 
E.Q, Ammann, M.K, Oshman, J.D. Foster, and J.M. Yarborough, "Repetitively 
Pumped Parametric Oscillator at 2.13 ,"" Appl. Phys. Letters 15, 131 (1969). 
E.0, Ammann and J. M, Yarborough, "Optical Parametric Oscillation in 
Proustite," Appl, Phys. Letters 17, 233 (1970). 
R.W. Wallace, "Stable, Efficient Optical Parametric Oscillators Pumped 
with Doubled Nd:YAG,"" Appl. Phys. Letters 17, 497 (1970). 
E£.0. Ammann, J.M. Yarborough, and J. Falk, "Simultaneous Optical Parametric 
Oscillation and Second Harmonic Generation," J. Appl. Phys. 42, S618 (1971). 
C.F. Dewey, Jr. and L.O. Hocker, "Infrared Difference-Frequency Generation 
Using a Tunable Dye Laser," Appl. Phys. Letters, 18, 58 (1971). 
J.M. Yarborcugh and G.A. Massey, "Efficient High-Gain Parametric Generation 
in ADP Continuously Tunable Across the Visible Spectrum," Appl. Phys. 


Letters 18, 438 (1971). 


4S 


a 


= tin er lll alka tl ci ee 


: 
i 


23. D.C. Hanna, R.C. Smith, and C.R, Stanley, "Generation of Tunable Mediun 
Infrared Radiation by Optical Mixing in Proustite," Opt, Commun, 4, 300 
(1971). 

24. D.C, Hanna, B, Luther-Davies, H.N. Rutt, and R.C, Smith, "Reliable 
Operation of a Proustite Parametric Oscillator,’ Appl. Phys, Letters 20, 
34 (1972). | 

25. L.S. Goldberg, "A Repetitively Pulsed Lil0, Internal Optical Parametric 
Oscillator,'' Presented at the VII International Quantum Electronics 
Conference, Montreal, Canada (1972). 

26. R.L. Herbst and R.L. Byer, "CdSe Infrared Parametric Oscillator,” 
Presented at the VII International Quantum Electronics Conference, 
Montreal, Canada (1972). 

27. R.L. Herbst and R.L. Byer, "Singly Resonant CdSe Infrared Parametric 
Oscillator,‘' Appl. Phys. Letters 21, 189 (1972). 

28. D.C, Hanna, B, Luther-Davies, and R.C, Smith, "Singly Resonant Proustite 
Parametric Oscillator Tuned from 1,22 ym to 8.5um," Appl. Phys. Letters 
22, 440 (1972). 

29, D.W, Meltzer and L.S, Goldberg, ‘'Tunable iR Difference Frequency Generation 
in Lil0,;" Opt. Commun. 5, 209 (1972). 

30. A.A. Davydov, L.A. Kulevskii, A.M. Prokhorov, A.D. Savel'ev, V.V, Smirnov, 


and A.V. Shirkov, "A Tunable Infrared F.rametric Oscillator in a CdSe 


Crystal,'' Opt. Commun. 9, 234 (1973). 
31. G.C, Bhar, D.C. Hanna, B. Luther-Davies, and R.C. Smith, "Tunable Lown- 


Conversion from an Optical Parametric Oscillator," Opt. Commn. 6, 523 


(1972). 


3 
q 
a | 
a 
. 4 
| 


i 32, P.P. Sorokin, J.J. Wynne, and J.R. Lankard, "Tunable Coherent IR Source 


Based Upon Four-Wave Parametric Conversion in Alkali Metal Vapors," Appl. 


hehe einstellen Oe 


b Phys. Letters 22, 342 (1973); Electro-Optical Systems Design 5, 4 (June 1973). 
46 


33. 


35. 


36. 


37. 


38, 


39. 


1 ETERNITY FAS 


40. 


41. 


42, 


43, 


nists she skis hn ania lerereeaas cies Ace Aa 2 te ren NTT A = 


C.D. Decker and F.K. Tittel, "High Power Broadly Tunable Difference 


Frequency Generation in Proustite," Appl. Phys. Letters 22, 411 (1973). 

D.C. Hanna, V.V. Rampal, and R.C. Smith, "Tunable Infrared Down Conversion 

in Silver Thiogallate," Opt. Commun. 8, 151 (1973). 

C.D. Decker and F.K, Tittel, "Difference Frequency Generation by Optical 
Mixing of Two Dye Lasers in Proustite," Opt, Commun. 8. 244 (1973). 

R.W. Wallace and S,E. Harris, Electronically Tunable Filter and Dye 

Laser, AFAL-TR-73-411, Nov 1973, 

R.L. Byer, M.M. Choy, R.L. Herbst, D.S. Chem’<, and R.S. Feigelson, 

"Second Harmonic Generation and Infrared M’.xing in AgGaSe.," Appl. Phys. 
Letters 24, 65 (1974). 

G.D. Boyd and D.A. Kieinman, "Parametric Interaction of Focused Gaussian 
Light Beams," J. Appl. Phys. 39. 3597 (1968). 

G.D. Boyd, H. Kaspar, and J.H. McFee, "Linear and Nonlinear Optical Properties 
of AgGaS., and CulnS, and Theory of the Wedge Technique for the Measurement 

of Nonlinear Coefficients," IEEE J. Quant. Electr. QE-7, 563 (1971). 

D.S. Chemla, P.J. Kupecek, D.S. Robertson, and R.C. Smith, "Silver Thiogallate, 
A New Material with Poter*ial for Infrared Devices," Opt. Commun. 3, 29 (1971). 
G.D. Boyd, H.M. Kaspar, J.H. McFee, and F.G. Storz, "Linear and Nonlinear 
Optical Properties of Some Ternary Selenides," IEEE J. Quant. Electr. QE-8, 
900 (1972). 

H, Kildal and J.C. Mikkelsen, "The Nonlinear Optical Coefficient, Phase- 
matching, and Optical Damage in the Chalcopyrite AgGaSe, ," Opt. Commun. 9, 

315 (1973). 


R.L. Byer, H. Kildal, and R.S. Feigelson, "CdGeAs.-A New Nonlinear Crystal 


2 
Phasematchable at 10.6 um," Appl. Phys. Letters 19, 237 (1971). 


47 


| ee Dee 


44, 


45. 


46. 


47. 


48, 


49. 


50. 


$1. 


52. 


53. 


54. 


G.D. Boyd, E. Beuhler, F.G. Storz, and J.H. Wernick, "Linear and Nonlinear 
Optical Properties of Ternary all plV cy Chalcopyrite Semiconductors," 
TEBE J. Quant. Electr. QE-8, 419 (1972). 

G.D. Boyd, E. Beuhler, and F.G. Storz, "Linear and Nonlinear Optical 


Properties of ZnGeP., and CdSe," Appl. Phys. Letters 18, 301 (1971). 


G.D. Boyd, T.J. Bridges, and E.G. Burkhardt, "Up-Conversion of 10.6 


Radiation to the Visible and Second Harmonic Generation in HgS," 

TEEE J. Quant. Electr. QE-4, 515 (1968). 

W.L. Bond, G.D. Boyd, and H.L. Carter, Jr., "Refractive Indices of HgS 
(Cinnabar) Between 0.62 and 11 u," J. Appl. Phys. 38, 4090 (1967). 

R.A. Andrews, "IR Image Parametric Up-Conversion," IEEE J. Quant. 

Electr. QE-6, 68 (1970). | 

D.S. Chemla, Ph.J. Kupecek, C,A, Schwartz, "Redetermination of the 
Nonlinear Optical Coefficients of Proustite by Comparison with Pyrargyrite 
and Gallium Arsenide," Opt. Comm. Vol. 7, 225 (1973). 

K.F, Hulme, 0, Jones, P.H. Davies, and M.V. Hobden, "Synthetic Proustite 
(Ag,AsS,): A New Crystal for Optical Mixing," Appl. Phys. Letters 10, 
133 (1967). | 

W.T. Gandrud, G.D. Boyd, J.H. McFee, and F.H. Wehmeier, "Ncnlinear 
Optical Properties of Ag,SbS,," Appl. Phys. Letters 16, 59 (1970). 

J.D. Feichtner, R. Johannes, and G.W. Roland, "Growth and Properties of 
Single Crystal Pyrargyrite," Appl. Optics 9, 1716 (1970). 

R.L. Byer and S,E. Harris, "Power and Bandwidth of Spontaneous Parametric 
Emission," Phys. Rev..168, 1064 (1968). . 

G.D. Boyd, W.L. Bond, and H.L, Carter, "Refractive Jndex as a Function of 


Temperature in LiNbO, ," J. Appl. Phys. 38, 1941 (1967) 


48 


ov aan ae Se TO URNS OO VCEEL ETS SEMITOSPUCUOUVOVOGET OTT OW OS RR oo 


tinea ieee as ct het wa ieMbhcnst te 5, 


1 

| 

4 

| 

\ 

q 

i 

: 
di 


r 


5S. J.D. Feichtner and G.W. Roland, "Optical Properties of a New Nonlineer 
Optical Material: T1,AsSe,," Appl. Optics 11, 995 (1972). 

56. A.J. Campillo and C.L. Tang, "Spontaneous Parametric Scattering of 
Light in Lil0,," Appl, Phys. Letters 16, 242 (1970). 

r> G, Nath and S, Hausstihl, "Large Nonlinear Optical Coefficient and 
Shase Matched Second Harmonic Generation in Lil0,," Appl. Phys. Letters 
14, 154 (1969). 

58. R.L. Herbst and R.L, Byer, "Efficient Parametric Mixing in CdSe," 
Appl. Phys. Letters 19, 527, (1971). 

59, J. Warner, "Phase-Matching for Optical Up-Conversion with Maximum 

a Angular Aperture - Theory and Practice," Opto-Electronics 1, 25 

(1969). 


49 


eee rmn OS aI 


| 
( 
| 


Mes ee ee 2 


APPENDIX A 
Ag, ASS, ( Proustite 


Sellmeier Bquation:“® 


tn°)? = 7.4822 + 0.4635 3 rebie 38 
r»* ~ 0.1160 
(n°)? 2 6.3434 + _ 0.3352 |g 007 32 
d? ~ 0.1117 
Refractive Index Data: 56,59 
A__(um) — (ua) n° n° 
0.5876 1.7018 2.7896 
0.6328 1.5803 3.0190 2.7391 
0.6516 1.5347 2.9970 2.7230 
0.6678 1.4975 2.9804 2.7094 
0.6943 1.4403 2.9550 2.6910 
1.014 0.9862 2.8264 2.5901 
1.129 0.8857 2.8067 2.8786 
1.367 0.7315 2.7833 2.5570 
1.530 0.6536 2.7728 2.5485 
1.709 0.5851 2.7654 2.5423 
2.50 0.4000 2.7478 2.5282 
3.56 0.2809 2.7379 2.5213 
4.62 0.2165 2.7318 2.5178 
10.6 0.0943 2.7030 2.5040 
50 
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-0.2799 
-0.2740 
-0.2710 
~0.2640 
-0.2363 
-0.2311 
-0.2263 
-0.2243 
-0.2231 
-0.2196 
-0.2166 
-0.2140 
-0.1990 
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AgGaS,, (Silver Thiogallate) 


Sellneier Bquation: “° 
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(n*)? «© . 48 + 0.2410 - 0.00210: a” 
d7 ~ 0.0870 
(n°)? = 5.497 + 0.2026 - 0.00283 2” 
2 ~ 0.1307 

Refractive Index Data:>” 

A_(ya) 7) (yan) n° n° n° - n° | 
fs | 
| .4900 2.0488 2.7148 2.7287 .0138 

.5000 2.9000 2.6916 2.6867 ~.0049 

.$250 1.9048 2.6503 2.6239 - 0264 | 

.5500 1.8182 2.6190 2.5834 -.0356 

.5750 1.7391 2.5944 2.5537 -.0407 | 

.6000 1.6667 2.5748 2.5303 -.0444 | 

.6250 1.6000 2.5577 2.5116 -.0461 

.6500 1.5385 2.5437 2.4961 -.0476 | 

.6750 1.4815 2.5310 2.4824 - 0486 

.7000 1.4286 2.5205 2.4706 -.0499 | 

. 7500 1.3333 2.5049 2.4540 -.0S09 

. 8000 1.2500 2.4909 2.4395 -.0514 | 

.8500 1.1765 2.4802 2.4279 ~.0522 

.9000 1.1121 2.4716 2.4192 -.0525 

.9500 1.0526 2.4644 2.4118 -.0526 | 

1.0000 1.0000 2.4582 2.4053 -.0529 i 
| 1.1000 .9091 2.4486 2.3954 - 0532 me 
1.2000 8333 2.4414 2.3881 -.0533 a 
| 1.3000 . 7692 2.4359 2.3819 -.0540 | 
. 1.4000 7143 2.4315 2.3781 -.0534 
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4 
eee) Satie ene Sal alae tm 
ee ee ere eee 2 


d (un) 


1.5000 
1.6000 
1.8000 
2.0000 
2.2000 
2.4000 
2.6000 
2.8009 
3.0000 
3.8000 
3.4000 
3.6000 
3.8000 
4.0000 
4.5000 
5.0000 
5.$000 
6.0000 
6.5000 
7.0000 
7.5000 
8.0000 
8.5000 
9.0000 
9.5000 
10.0000 
10.5000 
11.0000 
11.5000 
12.0000 
12.5000 
13.0000 


2.4280 
2.4252 
2.4206 
2.4164 
2.4142 
2.4119 
2.4102 
2.4094 
2.4080 
2.4068 
2.4062 
2.4046 
2.4024 
2.4024 
2.4003 
2.3955 
2.3938 
2.3908 
2.3874 
2.3827 
2.3787 
2.3757 
2.3699 
2.3663 
2.3606 
2.3548 
2.3486 
2.3417 
2.3329 
2.3266 
2.3177 
2.3076 
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2.3745 
2.3716 
2.3670 
2.3637 
2.3604 
2.3583 
2.3567 
2.3539 
2.3545 
2.3534 
2.3522 
2.3511 
2.3491 
2.3488 
2.3461 
2.3419 
2.3401 
2.3369 
2.3334 
2.3291 
2.3252 
2.3219 
2.3163 
2.3121 
2.3064 
2.3012 
2.2948 
2.2880 
2.2789 
2.2716 


-.0535 
- 0535 
-.0536 
-.9S27 
~ O287 
-.0535 
-.0535 
-.0535 
-.0535 
-.0534 
-.0840 
-.0535 
~.0533 
~.0536 
-,0542 
-. 0536 
- 0537 
v.0539 
-.0540 
- 0536 
~.0535 
-.0538 
-.0536 
~.0542 
-.0542 
-.0536 
-.0538 


-.0537 © 


7e 0540 
-.0550 


ST ESR AS eo: <a Seed Ser BS aS Powe Mae Vp a By Se ae 
Ree tee Cth ar ue 
“ a a | cae ma lh cha al EC NORA eneE d * ‘ee ws oot 


AgGaSe. 


Sellmeier Bquation : *? 
(n*)? = 3.9362 + 2.9113 ¥ 1.7954 
1- (0.38821 1-740 


(n) + ©=©6$.3332 + 3.3616 * 1.7677 


ca) 


i= (@.seaur$ 


e@ 
ww 
i) 
| 
> 
=] 
am ane -rearnserence ween Eta SIPS UTIeCAUSASANDIS ica SARE i EASUOEAESAEDNTE SeenSAOERY ct 


Refractive Index peta: 2 

d_(ua) a7) un”) n° n° n° - n° 

.7250 1.3793 2.8452 2.8952 .0480 

. 7500 1.3333 2.8191 2.8415 .0224 

. 8000 1.2500 2.7849 2.7866 .0017 

.8500 1.1765 2.7598 2.7522 -.0077 : 
.9000 1.1211 2.7406 2.7275 -.0130 

.9500 1.0526 2.7252 2.7085 -.0167 
1.0000 1.0000 2.7152 2.6934 -.0198 
1.1000 .9091 2.6942 2.6712 -.0230 
1.2000 .8333 2.6806 2.6554 -.0253 
1.3000 . 7692 2.6705 2.6438 -.0267 

1.4000 7143 2.6624 2.6347 -.0277 
1.6000 .6250 2.6516 2.6224 =v 0292 
1.8000 .5556 2.6432 2.6131 -.0301 
2.0000 .5000 2.6376 2.6071 ~.0305 
2.2000 .4545 2.6336 2.6027 -.0309 
2.4000 .4167 2.6304 2.5992 -.0313 1 | 
2.6000 . 3846 2.6286 2.5968 -.0317 } 
2.8000 3572 2.6261 2.5943 -.0318 
3.0000 . 3333 2.6245 2.5925 -.0320 
3.2000 .3125 2.6231 2.5912 -.0319 
3.4000 .2941 2.6221 2.5899 ~.0321 
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Aah nae 


b 
"4 
I 
; 
b 
‘ 
: 


A Cur) av? (um) 
5.6000 .2778 
$.8000 . 2632 
4.0000 . 2500 
4.5000 2222 
5.0000 . 2000 
5.5000 1818 
6.0000 . 1667 
6.S0C0 .1538 
7.0000 1429 
7.5000 1333 
. 8.0000 .1250 
8.5000 1176 
9.0000 lll 
9.5000 .1053 
10,000¢ . 1000 
10.5000 .0952 
11.0000 .0909 
11.5000 .0879 
12,0000 0835 
12.5000 -0800 
13.0000 0769 
13.5000 0741 


AAT TT Ae Rane OATRRINC REE Ae Ofori sae + - 


2.6213 
2.6200 
2.6189 
2.6166 
2.6144 
2.6128 
2.6113 
2.6094 
2.6070 
2.6049 
2.6032 
2.6009 
2.5988 
2.5964 
2.5939 
2.58917 
2.5890 
2.5868 
2.5837 
2.5805 
2.5771 
2.5731 
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2.5889 
2.5876 
2.5863 
2.5840 
2.5819 
2.5800 
2.5784 
2.5768 
2.5743 
2.5723 
2.5704 
2.5681 
2.5659 
2.5635 
2.5608 
2.5585 
2.5555 
2.5536 
2.5505 
2.5473 
2.5439 
2.5404 


ae Ms Sr 


eee doe el de 


_ ao POT 5 +) a erating i eg 
barter 


Sellmeier Equation: 


(n*)? 2 1 6 


Refractive Index Data: 


Ag,SbS, (Pyrargyrite) 


6.885 2 + ~—so0.ss a2 
a* . 0.16 r* . 295 


2 0.0202 x2 


=u A 
* - 0.16 x - 225 


Data not tabulated. 
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| 
! 
| 
| 


| 
| 


so aah wal 


MBAR iE Renn EER iain oi! eR nae ee 


Mt as ea i a ad et 


i 


(n®? = 4 


(n®)? = 4 


A (um) 


2. 3000 
2.4000 
2.5000 
2.6000 
2.7000 
2.8000 
2.9000 
3.0000 
3.1000 
3.2000 
3.4000 
; 3.6000 
i 3.8000 
4.0000 
4.2000 
4.4000 
4.6000 
4.8000 
5.0000 
5.5000 
6.0000 
6.5000 


Sellmeier Equation: 4° 
+ 8.89) 


Oey 


¢ 9.521 


+ Gey 


Refractive Index Data: “4 


av) (wavy 


- 4348 
+4167 
- 4000 
- 3846 
- 3704 
- 3571 
- 3448 
- 3333 
+ 3226 
~ 3125 
«2941 
- 2778 
- 2632 
- 2500 
- 2381 
+2273 
- 2174 
- 2083 


+ .2000 


- 1818 
- 1667 
- 1538 


n° 


3.6076 
3.5973 
3.5895 
3.5823 
3.5778 
3.5721 
3.5684 
3.5645 
3.5615 
3.5581 
3.5536 
3.5503 
3.5468 
3.5440 
3.5415 
3.5391 
3.5372 
3.5354 
3.5336 
3.5285 
3.5251 
3.5223 


1.909 


1 “( 


3.7545 
3.7316 
3.7186 
3.7030 
3.6926 
3.6846 
3.6775 
5.6714 
3.6661 
3.6574 
3.6508 
3.6454 
3.6402 
3.6368 
3.6329 
3.6299 
3.6273 
3.6249 
3.6178 
3.6134 
3.6104 


easement ae nl 


wR a te 


BE. MIRE OWS EP ea S ee eee. 


Sete See eea: © Sea ay 


7.0000 1429 3.$200 3.6078 08673 
7.8000 1333 3.$175 3.6050 0875 
8.0000 - 1280 3.5157 3.6030 0873 
8.5000 1176 3.5340 3.6009 0869 
9.0000 lll 3.5120 3.5988 0868 
9.5000 -1083 3.5098 3.5966 0867 
10.0000 . 1000 3.5078 3.5942 0864 
10.5000 0982 3.5054 3.5922 0868 
11.0000 - 0909 3.5031 3.5896 0865 
11.5000 0870 3.5004 3.5871 0868 
12.0000 0833 3.4977 
12.5000 0800 3.4980 
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(n*) ? 


A (ya) 
1.0139 
1.1287 
1.3673 
1.5395 
1.7109 
2.3253 
3.0 
4.0 
5.0 
6.0 
7.0 
8.0 
9.0 
10.0 
11.0 
12.0 


Sellmeier Equation: 


Refractive Index Data: 


31 


(n®)? © 4.1321 + 1.8587 A? 


A sal 0.2187 


© 4.0829 + 2.0038 27 


Av + 0.2075 


$8 


0.9863 
0.8860 
0.7314 
0.6538 
0.5845 
0.4301 
0.3333 
0.2500 
0.2000 
0.1667 
0.1429 
C.1250 
0.1111 
0.1000 
0.0909 
0.0833 


CdSe 
» 3.0461 22 
n*- 3380 
+ 3.5540 4 
v4. 3629 
n° n° 

2.5481 2.8677 
2.5246 2.5444 
2.4971 2.5170 
2.4862 2.5059 
2.4776 2.4974 
2.4627 2.4823 
2.4553 2.4748 
2.4500 2.4694 
2.4464 2.4657 
2.4434 2.4625 
2.4398 2.4586 
2.4367 2.4552 
2.4333 2.4534 
2.4294 2.4475 
2.4252 2.4430 
2.4204 2.4379 
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0,0196 
0.0198 
0.0199 
0.0198 
0.0198 
0.0196 
0.0195 
0.0194 
0.0193 
0.0191 
0.0188 
0.0185 
0.0181 
0.0181 
0.0178 
0.0175 


alla li, lan” Blt ae a 


i re eerie mm 


HgS (Cinnabar) 


Sellmeier Equation: “8 ; 
(n*)? = 6.9445 + _ 0.3688 5. Rae 
r=. 0.1387 
(n°) = 8.3922 + _ 0.5390 Bakunin? 
A” - 0.1388 
Refractive Index Date: 
A_(um) yt (att) n° ne n® - n° | 
i 
0.62 1.6129 3.2560 2.9028 0.3532 | 
0.65 1.5385 3.2064 2.8655 0.3409 3 
0.68 1.4706 3.1703 2.8384 0.3319 { 
0.70 1.4286 3.1489 2.8224 0.3265 | 
0.80 1.2500 3.0743 2.7704 0.3039 1 
0.90 1.1121 3.0340 2.7387 0.2957 d 
1.00 1.0000 3.0050 2.7120 0.2930 
1.20 0.8333 2.9680 2.6884 0.2796 
1.40 0.7143 2.9475 2.6730 0.2745 
1.60 0.6250 2.9344 2.6633 0.2711 
1.80 0.5556 2.9258 2.6567 0.2691 
2.00 0.5000 2.9194 2.6518 6.2676 
2.20 0.4545 2.9146 2.6483 0.2663 
2.40 0.4167 2.9108 2.6455 0.2653 
2.60 0.3846 2.9079 2.6433 0.2646 
2.80 0.3571 2.9052 2.6414 0.2638 
3.00 0.3333 2.9036 2.6401 0.2635 
3.20 0.3125 2.9017 2.6387 0.2630 
3.40 0.294" 2.9001 2.6375 0.2626 
3.60 0.2778 2.8987 2.6358 0.2629 
3.80 0.2632 2.8971 2.6353 0.2618 
4.00 0.2500 2.8963 2.6348 0.2615 
5.00 0.2000 2.8863 2.6267 0.2596 
6.00 0.1667 2.8799 2.6233 9.2566 
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i ok ee ee in eh teehee oe te ek i 


7.00 
8.00 
9.00 
10.00 
11.00 


0.1429 


0.1250 
9.1111 
9.1000 
0.0909 
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Semmens oT 


'2.6156 


2.6112 
2.6066 
2.6016 
2.5914 


e 
nm -n 
ae 


0.2585 
0.2562 
0.2542 
0.2504 
0.2520 


Fhe al Raat ita ais I taal, 


: a a Le eC eT ee TOUTE 
Mee ate et tein, a ee ae a 


POD Se hm hls Bk oa te 


. 
hen. tee 


Sellmeier Bquation:* 


(n*)? = 3.463 + 0.0224 - 0.00677 27 
A* ~ 0.0958 
6.2 2 
A” ~ 0.0915 
57 
Refractive Index Data:°°* 
- a | 
A_(um) A 1 (um ) n° n° n° - n° 
0.436 2.294 1.931 1.766 -0.165 
0.500 2.000 1.908 1.754 ~0.1S54 
0.530 1.887 1.901 1.750 -0,151 
0.578 1.730 1.888 1.742 -0.146 
0.690 1.449 1.875 1.731 -0.144 
0.800 1.250 1.868 1.724 -0.144 
1.060 0.943 1.860 1.719 -0.141 


*Equation fit to refractive index data anc curves in Ref. S6 and 57, 
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LiNbO, 
| Sellmeier Equation: “® 
| (n°)? = 4.9260 + 0.1170 - 0.0275 a2 
7 . 0.0473 
(n°)? = 4.5778 + 0.0964 - 0.0221 4? 
r - 0.0446 
Refractive Index Data:°* 
A_(ua) a7) um) n° n° n° - n° 
0.42 2.3810 2.3038 2.4144 -0.1106 
: 0.45 2.2222 2.2765 2.3814 -0.1049 
0.50 2.0000 2.2446 2.3444 -0.0998 
0.55 1.8182 2.2241 2.3188 -0.0947 7 
0.60 1.6667 2.2083 2.3002 -0.0919 
0.65 1.5385 2.1964 2.2862 - 0.0898 
0.70 1.4286 2.1874 2.2756 - 0.0882 
0.80 1.2500 2.1741 2.2598 -0.0857 
0.90 1.1111 2.1647 2.2487 -0.0840 
1.00 1.0000 2.1580 2.2407 * 20.0827 
1,20 0.8333 2.1481 2.2291 -0.0810 
1.40 0.7143 2.1410 2.2208 -0.0798 
1.60 0.6250 2.1351 2.2139 -0.0788 
1.80 0.5556 2.1297 2.2074 -0.0777 
2.00 0.5000 2.1244 2.2015 -0.0771 
2.20 0.4545 2.1187 2.1948 -0.0761 
2.40 0.4167 2.1138 2.1882 -0.0744 
2.60 0. 3846 2.1080 2.1814 -0.0734 
2.80 0.3571 2.1020 2.1741 0.0721 
3.00 0.3333 2.0955 2.1663 -0.0708 i 
3.20 0.3125 2.0886 2.1580 -0.0694 : 
3.40 0.2941 2.0814 2.1493 -0.0679 : 
3.60 0.2778 2.0735 2.1398 0.0663 
3.80 0.2632 2.0652 2.1299 . -0.0647 
4.00 0.2500 2.0564 2.1193 -0.0629 | 
62 4 


T1 ,AsSe, 


Sellmeier Equation: > 


(rn)? © 1 + 9.97707 + 0.067 2" 


12. (0.435) rA“e (20) 


(ny? 2 1 + 8.78507 + 0.051 2 


d2. (0.438)* “= (20) 


Refractive Index Data o 


A_ (uit) acl qua”) n° n° 

1,583 t 0.03 0.644 3.443 3.248 
2.66 + 0.05 0.376 3.356 3.170 
3.29 +0.1 0.304 3.339 3.152 
3.365 + 0.065 0.297 3.337 3.155 
3.38 + 0v16 0.296 3.339 3.152 
4.35 + 0.09 0.230 3.332 3.148 
4.46 +£0.17 0.224 3.33" 3.142 
4.55 + 0.2 0.220 3.326 3.142 
§.26 +0.3 0.190 3.321 3.141 
5.3 2£0.1 0.189 3.326 3.142 
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-0.198 
-0.186 
-0.187 
-0.182 
-0.187 
-0.184 
-0.192 
-0.184 
-0.180 
-0.184 


sand Salih ie 23. ai 


Sellmeier Equation: - 


(n*)? = 9.7497 + 0.6889 


A” - 0.1402 
(n°)? = 9.9870 + 0.7451 

A“ - 0.1489 
Refractive Index pata: 45 
A_(ua) a7! (um) n° 
0.64 — 1.5625 3.5052 
0.66 1.5152 3.4756 
0.68 1.4706 3.4477 
0.70 1.4286 3.4233 
0.75 1.3333 3.3730 
0.80 1.2500 3.3357 
0.85 1.1765 3.3063 
0.90 1.1111 3.2830 
0.95 1.0526 3.2638 
1.00 1.0000 3.2478 
1.10 0.9091 3.2232 
1.20 0.8333 3.2054 
1.30 0. 7692 3.1924 
1.40. 0.7143 3.1820 
1.60 0.6250 3.1666 
1.80 0.5556 3.1562 
2.00 0.5000 3.1490 
2.20 0.4545 3.1433 
2.40 0.4167 3.1388 
2.60 0.3846 3.1357 
2.80 0.3571 3.1327 
3.00 0.3333 3.1304 
3.20 0.3125 3.1284 
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0.00278 27 


0.00271 a° 


3.5802 
3.5467 
3.5160 
3.4885 
3.4324 
3.3915 
3.3593 
3.3336 

3.3124 
3.2954 
3.2688 
3.2493 
3.2346 
3.2244 
3.2077 
3.1965 
3.1889 
3.1829 
3.1780 
3.1745 
3.1717 
3.1693 
3.1671 


6 
n -?R 


0.0750 
0.0710 
0.0684 
0.0652 
0.0595 
0.0558 
0.0530 
0.0596 
0.0486 


0.0476. 


0.0456 
0.0438 
0.0423 
0.0423 
0.0411 
0.08603 
0.0399 
0.0396 
0.0391 
0.0388 
0.0390 
0.0388 
0.0386 


ve ap ee OE CM. TS OD ee 


~~. a lime ES ttt 


wc 


3.40 
3.60 
3.80 
4.00 
4.20 
4.50 
4.70 
5.00 
5.50 
6.00 
6.50 
7.00 
7.50 
8.00 
8.50 
9.00 
9.50 
10.00 
10.50 
11.00 
11.50 
12. 00 


7) ie) 


0.2941 
0.2778 
0.2632 
0.2500 
0.2381 
0.2222 
0.2128 
0.2000 
0.1818 
0.1667 
0.1538 
0.1429 
0.1533 
0.1250 
0.1176 
0.1111 
0.1053 
0.1000 
0.0952 
0.0909 
0.0870 
0.0833 


n° 


3.1263 
3.1287 
3.1237 
3.1223 
3.1209 
3.1186 
3.1174 
3.1149 
3.1131 
3.1101 
3.1057 
3.1040 
3.0994 
3.0961 
3.0919 
3.0880 
3.0836 
3.0788 
3.0738 
3.0689 
3.0623 
3.0552 


3.1647 
3.1632 
3.1616 
3.1608 
3.1595 
3.1561 
3.1549 
3.1533 
3.1518 
3.1480 
3.1445 
3.1420 
3.1378 
3.1350 
3.1311 
3.1272 
3.1231 
3.1183 
3.1137 
3.1087 
3.1008 
3.0949 


*Equation fit to refractive index data in Ref. 45. 


0.0384 
0.0376 
0.0380 
0.0386 
0.0386 
0.0374 
0.0375 
0.0383 
0.0387 
0.0379 
0.0387 
0.0380 
0.0384 
0.0389 
0.0392 
0.0392 
0.0395 
0.0396 
0.0399 
0.0398 
0.03586 
0.0397 
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APPENDIX B 


The calculations presented in this report were carried out on a CDC 


6600 computer, using a program called NONLIN. A listing of program NONLIN 
is given in this Appendix, 


Three units must be added to NONLIN in order to perform a calculation. 
The first unit consists of the index of refraction equations for the material 
of interest. The indices are defined by the arithmetic function OINDEX (X) 
and EINDEX (X), when X is wavelength in microns. These two functions are 
to be inserted before the first executable statement in NONLIN; a pair of 
comment statements (see lines NON 74 and NON 75) have been included in the 
program to indicate the correct location for the index functions. (In the 
actual deck of cards used, this pair of ~omment cards was of a different 
color than the rest of the deck to aid in locating them.) A listing of the 
index of refraction functions for all the materials appearing in Table III 
is given in Appendix C. 


hain nino. ee ee neo ee 


The second unit that must be added to NONLIN is the subroutine DSQUAR. 
This subroutine is to be inserted immediately after NONLIN. The purpose of 
DSQUAR, as oxplained in the comments section of NONLIN, is to compute the 
square of the effective nonlinear optical coefficient as a function of the 
phase matching angle ©. This subroutine will vary depending on the crystal- 
lographic point group to which the material of interest belongs. The various 
DSQUAR subroutines that were used for this report are listed in Appendix D. 


These subroutines are based on the expressions for dire (9) given ‘n Table III, 


The third unit to be added to NONLIN is the input. For the CDC SCOPE 


System, an end-of-record card (7/8/9 in column 1) must be inserted between 
subroutine DSQUAR and the input. 


As an example of the calculations done for .his report, the output 
: generated for cinnabar (HgS, point group 32) is reproduced in Appendix E. 
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HAOOMSPAIVOIOIOAIAPIGAGAOARASGAAGOGOAHAASHGOSHAHAAAGAGAHAHAIAIAIAGAGAGAIASAHAIAAARHAADOHHOS 


PROGRAM NONLIN( INPUT OUTPUT , TAPE SUINPUT, TAPERSOUTPUT) 


PROGRAM TO COMPUTE PHASE NATCHING CURVES AND SMALL SIGNAL GAIN FOR 
COLXNEAR THREE WAVE OPTICAL PARAMETRIC INTERACTION. THEORY OF 
CONFOCALLY FOCUSEN PHASE MATCHING IS USED (REFERENCES - ~ G.0. 80Y0 
AND 0. A. KLEINMAN (FR ANT! KD), “PARAMETRIC INTERACTION OF FOCUSED 
GAUSSTAN LIGHT BEANS™, J. APPL. PHYS. VOL. 3943597(1968), ANDO S$. Eo 
HARRIS. “TUNABLE OPTICAL PARAMETRIC OSCILLATORS”, IEEE PROC. VOL. S?y 
2096(1969)). A CRYSTAL THICKNESS OF 1 CH IS ASSUMED. 


TNPUT = = CARD 1 HEAN FORMAT(8A10) 
HEAD TS ANY 80 COLUMN COMMENT, 
WHICH WILL BE LISTED WITH THE 
NuTPur, 
CARD 2 THIF.O FORMAT(I3,5F 10.3) 


TECIAIF.LE.0) THE MATERIAL TO AE PHASE 
MATCHED IS NEGATIVE RIREF. 

TFCIBIF.GT.09 THE MATERIAL IS POSITIVE 
BIREFRINGENT. 


) ARE THE ABSOLUTE VALUES OF THE 
NLO COEFFICXENTS IN UNITS OF 
10*712 M/V, 


CAROS SebeFTC. PUMP,CUTOFF FORMAT(2F10.3) 


PuNP TS THE PUMP WAVELENGTH (CIN 
MICRONS). 
CUTOFF TS THE INFRARED LONG WAVE= 


LENGTH CUTOFF FOR THE TOLER, 
IN MICRONS. IT IS NEEDED TO 
PREVENT CALCULATION OF 
REFRACTIVE INDICES BEYOND THE 
RANGE OF VALIOITY OF THE INOEX 
EQUATIONS. AS MANY PUMP 
WAVELENGTHS AS DESIRED MAY BE 
TNCLUDED. ONE PER CARD. 


YHE OUTPUT IS SELF EXPLANATORY. THETA ITS THE P.M. ANGLE, RHO IS THE 
BIREFZINGENCE WALKOFF ANGLE. AND O-EFF IS THE EFFECTIVE NLO 
COEFFICIENT (VARTES WITH THETA). HAD ITS THE FUNCTION DEFINED AY A 
AND Ky (EQ 3,49). THIS FUNCTION TAKES WALKOFF INTO ACCOUNT THROUGH THE 
PARAMETER 3, OEFINED SY A AND K (EQ 3.35). GAIN/WATT (MAX) IS THE 
GAIN/WATT FOR AN INFINITELY THICKse LOSSLESS CRYSTAL. AREA IS THE 
OPTIMUM REAM WATST CROSS SECTION. 


THE SUBROUTINE NSQUAR(DSQ,THETAR.T) MUST 9E SUPPLIED AY THE USER. IT 
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WON 
WON 
WON 
WOR 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
WON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 


Ber ARF ane 


7 OSRAPHETSS “TERRES FF 


ot AE, 


> stn “ree etna. 


ada LO kek 


HOHONM AAA aooogdca 


RETUPNS THE VALUE OF D-EFF (SQUARED) AS A FUNCTION OF THETA, FOR THE 
FOLLUWING CASES - = 


T2Q NEGATIVE SYIREF., TYPE TIT P.M. (0 + O & E) 
T#i NEGATIVE AIREF., TYPE IT P.M. (9 + E 2 ED OR (FE #0 ® E) 
Te2 POSITIVE AIREF., TYPE T P.M (€ + © & 0) 
T#3 POSITIVE AIREF., TYPE IT P.M. (0 ¢ E = OP OR (FE + 0 ® 0) 


SUBROUTINE DSQUAR FUST ALSO CONTAIN THE FOLLOWING CARD 
COMMON/COEF/N(5) 

INSERT BETWEEN THE TWO BLUE COMMENT CARDS BELOW THE FUNCTIONS 

OINDEX(X) AND EINNEX(X), WHERE OINDEX AND EINDEX ARE FUNCTIONS THAT 

CALCULATE THE ORDINARY AND EXTRAORDINARY REFRACTIVE INDICES FOR THE 

MATERTAL AS A FUNCTION OF WAVELENGTH (X = WAVELENGTH IN MICRONS). 


OIMENSION HEADS) 
REAL TOLER LHS 
DATA NEGe IPOS eS 2aC3/BHNEGATIVE o BHPOSI TIVE +1. 189815528E-3, 
13.1415926°536/ 
DATA TAsTSeARgABG ACAD AFGAF/1H etHTetQHO ¢ O = E),10HO + E = Edy 
Z1JHE * O = EV VLOHE * E = 0)e10HO + E = OVF10HE + 0 = 00/ 
INSERT OINDEX(X) AND EINDEX(X) HERE. 
NOW YOU ARE IN BUSINESS. 
€12160./C3 
REAN(S,10) HEAN,IBIF.D 
10 FORMAT (BA10/T30 5F 10. 3) 
IF CISIF) 215515016 
15 IKINO®NEG 
GO TO 20 
16 IKIND=TPOS 
20 READ(3,25) PUMP CUTOFF 
25 FORMAT(2F 10.3) 
IF CEOF (39) 26427 
26 STOP 
27 PRINT 30 eHEAN eI KINO» PUNP oD, CUTOFF 
30 FIRMAT(1Hie*PROGRAM NONLIN=--COLINEARLY PHASE MATCHED OPO*/1HO/1H « 
18AL0/1HO/1H o*PHASE MATCHING FOR *%,A8,4* BIREFRINGENT CRYSTAL*/1HO, 
2°9UMP WAVELENGTH (MICRONS) «® oF12.3/1H0 oe NONLINEAR COEFFICYTENTS (X10 
3E1L12 M/VI EF, SFIS .S/IHOe*CUTOFF WAVELENGTH=*,F12.3) 
CUTOFF 81. /CUTGFF 
L=0 
T=0 
OMEGAP21./PUMP 
QEL =OMEGAP/150. 
XOZO INDEX (PUMP) 
31 SIGNAL22.*PUMP 
OMEGAS*16/SIGNAL 
OMEGAT=OMEGAS 
TOLERSSIGNAL 
IF(L.EQ.0) 60 TO 34 
GO TO (112611221060 41046) ot 
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NON 
NON 
NON 
NON 
NON 
WON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 
NON 


99 


NON100 
NON1O1 
NON102 
NON103 


4 


Sea eee oe 


3& TFECIBIF) 35.35,1000 
38 KESE INDEX (PUMP) 
YOuOINDEX (SIGNAL) 
77ERN2YO 
¥2XKO*FXKO-KESXKE 
TECK) Foe? Ooh 
eO IFCKE.LE.YO) GO TO 65 
PRINT 465.TA AA, PUMP 
&S FORMATC(1HOe*THIS MATERTAL CANNOT RE TYPE [T*eA1e2H (eAL0e% PHASE MA 
4TCHEN FoR THE 2UMP WAVELENGTH? .F1Z.3e% MICRON*) 
S30 OMEGASZOMEGASONEL 
IF (OMEGAS.GEeOMEGAP) GO TO 100 
OMEGATZCNEGAP-OMEGAS 
TF COMEGAY.«LTeSUTOFF) GO TO 100 
SIGNAL2=1./0MEGAS 
TNLERs1 . /OMEGAT 
TFC ITeEQe2) GO TO 1005 
YIZOINOEX (SIGNAL) 
PFOROINDEXCYNLER) 
TFC COMEGAP*YXED oGT. (COMEGAS*YOCOMEGAT*®720)) GO TO 56 
PRINT GO0eSIGNAL», IALER 
6G FORMAT(IH «*IN THE WAVELENGTH REGION*.F12.3¢* TO*eF12e3—% MICRON*) 
G9 TO 80 
Te PIINT 7SeITKINNeSTGNALSTOLER 
TS FORMATCLHOe* THE CALCULATED INDICES FOR THIS MATERIAL 0O NOT GIVE * 
1e4Be* ATREFRINGENCE FOR SIGNALEFeF1 203% sIOLER=*eF 12.34% NEXT CAS 
2€*) 
GO TO 110 
6§ 79=Y0 
80 PRINT 92eTA,AA 
Oi FIRMATCIHIe*RESULTS FOR TYPE T%eAle2H FeAL0e*% PHASE MATCHING ARE*/ 


11H) 
PRINT 90 

90 FORMATCIHOe*SIGNAL INMLER THETA RHO THETA RH 
10 O-EFF 8 H(A) GAIN/WATT GAIN/WATT ARE 


2A*/1H 11K o2CEK eSHIRADD Dd o 3K o2CIX SHINES) D 2% oH (FLO0EL2) o2HX,y 
BSHC ICM) BX eSHIMAYD 4 7Xe7H(SQ CM) /1H ) 
BS Y=(OMEGAP*XOZ (OMFGAS*YO ¢OMEGAT* 70)) P¥ 2-16 
SIN THEXE*SART CY/X%) 
IFCSINTH.GE.te) GO TO 1035 
THE TARZASTNISINTH) 
SINSQZSINTH*SINTH 
Ya *SINSO+KEFXE 
RHORZ.5*SIN( 2 o* THETAR) *X/V 
69 72XO*XE/SART(Y) 
88 RHOR=ATANCRHOR) 
THE TADSC1* THE TAR 
RHON=C1*RHOP 
CALL ODSQUARINS I .THETARST) 
R2Z7%(106.%RHOR) ¥*2/P UMP 
H121.068/ (1.41.068*9) 
H421./77 
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NONTO4 
NON105 
NON106 
NON107 
NON1LO8 
NON109 
WONT10 
WONL11 
NON112 
WON113 
NON114 
NONL15 
NON116 
NONL17 
NONL18 
WON119 
NON120 
NON121 
NON122 
NONL23 
NONI2& 
NON125 
WON1L26 
NON 27 
NONL28 
NONL29 
NON1L30 
NON131 
NONL32 
NON1LS3 
NONL34& 
NONL3SS 
NONLS6 
NON1L37 
NON138 
NON1L39 
NON140 
NON1L41 
NONL&2 
NON143 
WON1L4S 
NONLSS 
NONL46 
NON147 
NON14&A 
NONI49 
NON1SO 
NONLS1 
NONL52 
NON153 
NON1LS4& 
NON1I55 


RE UFR reper gay 


RAzC3*A/4, WON1S6 

IF(8.LE.2.) GO TO 93 NONLS?7 

R)eA/2, NONLS6 

69 TO 94 WONLS9 

93 Riz. WONL60 

9& AeSQRT(Aa) NONL62 

AREAZS .£-5#A029UMO/7 ZERO NONL62 

DEFFRSQRTINSQ) NON1L63 

GAINL2C240SQ*PUMP/ (STGNAL* SIGNAL *IOLER* TOLER *7*% 72) NONL6& 

GAINMZGAIN1 FHM NONL6S 

GAINIZGAINIFHL NON1L66 

PRINT BSeSIGNAL eS OLER, THETAR, RHOR, THE TAD eRHON NEFF o Betd eGAINIs NONL6? 

IGAINMSAREA NONL68 

95 FORMATULH FS e2eF Be 2elPELce SeEL Se 3eOPFI02e FBe2e F920 FI 2eF1Le Se NONI69 

13(21PE13.3)) NONL7O 

96 OYEGASZSOMEGASF9OEL NON172 

TF COMEGAS.GE.OMEGAP) GN TO 110 NONL72 

OME GATSOMEGAP ~OMFGAS NONL73 

IE (OMEGALT eLTe CUTOFF) GO TO 110 NONL74 

SIGNAL21./0MEGAS NONL7S 

¢. TILER=S1./9MEGAT NON 7A 

; VIZOINOEX (SIGNAL) NONL7? 

FIZOINDEXCTALER) NONL78 

YF(I.EQ.6) GO TO A5 . WONLT9 

GY TO (C108 01801650 oe LOMO LISI oL NONI 80 

100 PRINT 108 - NONLB4 

105 FIRMAT(iH ,* TYPE It TS ALSO EXCLUNED.*) NON1 82 

G) TO 20 NON183 

106 7©=sEINDEX( TOLER) NONLS4 

69 To 135 NONI85 

110 YFC(L.EQ.2) GO TO 20 NONL86 

TFCI.EQ.CY GO TO 11t NON18? 

TE(ILEQ.1) GO TO 160 NONL88 

TFE(L=6) 1955.108%,20 NON189 

142 L212 NONIS0 

T21 NONI92 

G? TO 31 NON192 

412 VOZOINDEX(SYIFNAL) NON193 

792Y0 NONL94 

FEREINDEYC(INLER) NONL9S 

TF CCYOO7E) GEol2cFXEV) GO THCL3IO01 75) ol NON196 

GD TO (1149164) of NON197 

416 PRINT 45,T%.AR.PUMO NON198 
: 12! NMEGASZOMEGAS 4D EL NON199° i 
TF COMEGAS.GFE2OMEGAP? 6D TO 160 NON200 ' 
OME GAIT ZOMEGAP=-OMFGAS NON202 { 
TF COMEGAT«LTeCUTNFF) GO TO 160 * NON2Z02 . 
SIGNALT1./0MEGAS NON203 \ 
VOSOINDEX (SIGNAL) NON264 
TNLER=1./9MEGAT NON?2CS 
: JEZEINDEX(IANLER) NON206 ; 
: IF ( COMEGAS* YO*OMFGAT®7E) LT. (CONEGA2*XE)) GO TO 120 NOH267 
i | 
70 ; 


130 
138 
137 
140 
ivi 
142 
145 


150 


155 


156 


169 


164 
176 


175 
180 


190 
195 


EEE Sea RT Ree RS ae a 


PRINT 60 eSIGNALAINLER 
ZISOTNDEX(TNLER) 

PRINT 91.TA,AR 

PRINT 90 

¥IZ70*Z0-7E*7E 

TE CNID 7H_7%hei37 

SINS4Q20. 

NEGse.i 
LHSROMEGAS*VOCOMEGAL ®70*®7E/SART (XT#SINSOF7ZE® Z7E) 
RHSBOMEGAP# XYOFXESSARTL MESINSQSVEF KE) 
TF(LHS.GE.°HS) GO TO 150 
SINSOSSINSONEC 

IF CSTNSQeGEeied GO TO 1035 

GO TO (140519501025) aL 

TF (DEC. LE et eF 6) GO TO 155 
STNSQZSSINSQ-NEC 

NEC#.1°0EC 

69 TO 146 

THE TAR2SQRTISINSO) 

THE TARZASIN(THE TAR) 
SINZ2.5*SIN(2.* THETAR) 
TE(L.EQ.3) GA TO 1045 
Y=X*SINSQ¢KEPXE 

RPSSIN2*X/Y 

GO TO (1564200) eL 
VYIT=EXIFSINSQA47E% 7E 
RIasSINe*yIsv_t 

PHORZAMAX!S (R2P.RPT) 

Gd Ta 89 

L22 

G93 TO 31 

PRINT 452TH ACePLIMP 
OVEGASZOMEGAS CNEL 

TF (COMEGAS.GE.OMEGAP) 69 TO 20 
OME GATS OMEGAP-OMFGAS 

JF COMEGAT.LT CUTOFF) GO TO 20 
SIGNAL21./70MEGAS 
VEZEINDEX(STGNAL) 
INLERs1./0NEGAT 

7FOZOINDEY (TOLER) 

TF CCOMEGAS* VE *CMFEGAT #70) LT. (OMEGAP*XE)) GO TO 179 
PRINT BO eSIGNAL eINLER 
YOZOINDEX (SIGNAL) 

PRINT 914-T9,4f 

PRINT 90 

YESEINOEX (SIGNAL) 
X¥SZYO*®VYO-VE*VE 

IF (XS) 74.746,190 

SINSQ20. 

OECzel 

LHS SOMEGAS*YOFVE/SIRTCKSFSINSOSVEFVE) COMEGAT F70 
G9 TO 141 


71 


Se ee ge ee 


menial 


ew Seer Ome 


NON20 8 
NON209 
NON210 
NON221 
WON212 
NON213 
NONR1LS 
NON215 
NON226 
NON217 
WON218 
NON219 
MON226 
NON221 
NON222 
NON223 
NON224& 
NON225 
NON226 
NON227 
NON228 
NON229 
NON230 
NON231 
NON232 
NON233 
NON234& 
NON2 35 
NON236 
NON237 
NON238 
NON239 
NON240 
NON242 
NON2Z42 
NOK243 
NON2Z&4& 
NON245 
NON246 
NON247 
NON248 
NON24&9 
NON250 
NON251 
NON252 
NON253 
NON254 
NON2S5 
NON2SE 
NON257 
NON2Z5S 
NON259 


1 
q 
4 
i 
{ 
{ 
4 


+ et gma 


1000 


1008 


1010 
1015 


1016 
1620 
1022 
1025 


1035 
1040 


1065 


1050 


1055 
1060 


YS2XS*SINSOSYVES VE 
RS#SIN2*xS/vs 
RHORRAMAXL(RE.RS) 

69 TO 89 

Y22 

23 

VESEINDEX(SIGNAL) 
TZERORVE 

TFCXO.LE.VE) GO TO 1020 
PRINT &5,TAcAN, PUMP 

GO To 50 

VEBEINDEX (SIGNAL) 
JFESEINUEXC(IALER) 

IF (C COMEGAP*# X0) GT. (OMEGAS®VESONEGAT*7E)) GO TO 80 
PRINT BOeSIGNAL TOLER 
69 TO 1018 

7TEZVE 

VOSOINOEX(SIGNAL? 
ZIOZOINDEXCTOLER) 

PRINT 91.TA,AN 

PRINT 90 

XSZVE*VYE-YO*YO 

TFE(XS) 76474-1020 
XTBZE*Z7E-70*70 

IF CKI) Fhe PHeiD2t 
SINSQs0. 

NECe#.4 

LHSZOMEGAS* YN*VE/SORT(IVE* VE -XS*SINSQ) SOMEGAL*®2Z0*7E/SORT (ZE% ZE= 


1 XT*SINSO) 


RYSBOMEGASD FXO 

GO TO 142 

PRINT 1040 ¢STGNAL,INLER 
FORMAT(1HO,*PHASE MATCHING CUTS OFF AT STGNAL®* ,F22.34% AND IOLER= 
1*.F12.3) 

G9 TO 96 
YVS2VE*FVE=xXS*SINSA 
VISZE*7E-XTFSINGO 
PS=SIN2*XS/YS 
RISSIN2*XT/YT 
PHORZAMAX1 (RS RT) 

7=Xx0 

60 TO 88 
YESEINOEX(STGNAL) 
7EZEINDEX(IALER) 

GO TO 1016 

T=3 

69 TO 31 

L=4 

YIZOINDEX (STGNAL) 
7EZEINDEXCIOLER) 

79270 

TFC (CYO*ZE) GE e(2.*X0)) GO TO 10790 
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WON260 
NON261 
NON262 
NON263 
NON2Z66 
NON265 
NON266 
NON267 
NON268 
NON269 
NON270 
NON271 
NON272 
NON273 
NON276 
NON275 
NON276 
NON277 
NON27 4A 
NON279 
NON280 
NON261 
NON282 
WON283 
NON2 84 
NON285S 
NON2Z86 
NON2Z47 
NON268 
NON2Z89 
NON290 
NON291 
NON292 
WON293 
NON296 
NON295 
NON296 
NON?297 
NON298 
NON2 99 
NON300 
NON301 
NONS02 
NON3O3 
NONS 04% 
NONSOS 
NON306 
NON307 
NCN308 
NON3SO9 
NON310 
NON311 


PUTER PRN Cysts Pat aia oe See ee Se 
or a ea ; 


1 2 ANN eRe 
PRINT 656 TRAE, PUMP NONS12 
1065 OYEGASEOMEGAS4NEL NONS23 
TFCOMEGAS.GE.OMEGAP) GO TO 31 NONS16 
OWEGATZOMEGAP-ONEGAS NONS15 
TF(OMEGAT.LT.CUTOFF) GO TO 31 NONS16 
SIGNAL#14/0MEGAS NONSLT 
TNLERe1  SOMEGAT NONS18 
YISOINOE X (SIGNAL) NONS29 
7EREINOEX (TOLER) NON320 
TFC COMEGAS* YNCOMEGAT*7E) «LT. (OMEGAP*XO)) GO TO 1065 NON324 
7IBOINOEX (TOLER) NON322 
PRINT 60eSIGNAL yTNLER NONS23 
1070 PRINT 91,7, AE NONS26 
PRINT 90 NON3S25 
1075 XTeZ£e7E-7NF70 NON326 
TE(XT) 74574610 76 NON327 
1076 V21.-(OMEGAT®70/ (OME GAP®X0-OMEGAS*YO) ) #42 NON328 
SINTHAZE*SORTUY/¥T) NON329 
IF(SINTH.GEs1.) GO TO 1035 NON330 
THE TARZASIN(SINTH) NON334 
SINSQ2SINTH*SIN TH NONS32 
VRZE*2E-XT#SINSO NONS33 
RHORS.S*SIN(2,* THETAR) *XISY NON3S34 
7=X0 NON335 
Go TO 88 NONS36 
1080 7EZEINDEX (TOLER) NON337 
60 YO 1075 NONS38 
1085 GO TO 34 NONS39 
1086 L=5 NONS40 
YESEINOEX (STGNAL) NIN341 
7IROINDEXCIOLER) NONSG2 
¥o270 NONS43 
TF CCVE¢70) GE. (2*X0)) GO TO 1095 NON344 
PRINT &5,TRAF, PUMP NON345 
1090 OVEGASZOMEGASENEL NONS46 
IF (OMEGAS.GE.OMEGAP) GO TO 20 NON347 
OMEGAT2OMEGAP-OMEGAS NONS4S 
IF (OMEGAT.LTsCUTOFF) GO TO 20 NON349 
SIGNAL: 1./0MEGAS NON3S0 
TNLER21./0MEGAT NON351 
YE®E INDEX (SIGNAL) NON352 
702 0IXNDEX (TOLER) NON353 
IF ( COMEGAS*YE+OMFGAT*70) LT. (OMEGAP*X0)) GO TO 1090 NONS54 
YO2OINDEX (SIGNAL) NON3SS 
PRINT 60 «eSTGNAL «TOLER NON356 
1095 PRINT 91.19,AF NON357 
PRINT 90 NON3S8— 
1096 XStYE*VE-YO*YO NON359— : 
TF (XS) The7GoA1Cd NONS60 
1100 Y=ie-(OMEGAS*YO/ (OMFGAP*KO-OMEGAT®70) ) #*2 NON361 | 
SINTH2YE*SORTIV/¥S) NON362 
TF(SINTH.GT.2.) GO TO 1035 NON363 
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THETAR=ASIN(SINTH) NON364 
| STHSQZSINTH*SINTH NON365 
VaYE#YE=XS*SINSO NON?66 
RHOR=.5*STN (20% THETAR) #XS/Y NON367 

72x0 NON368- 

G) TO 88 NON369 

1105 VEZEINDEX(SIGNAL) NON370 

GO TO 1096 NON374 
: END NON372 


j 
; 
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APPENDIX C 


TONE SS ake apne eee rae: TY Fare 


er ee Se ae ott 


PETS PE Se ae eee et eee ee aS Me RO TAM See ge eee See tere CAEL Gee 


SILVER GALLTUM SULFINE 
OTNDEX(K)=SOOT(E, 7284.241 SIK*X=.087 1-.0021 FX*X) 
FINDEX(X)=SART (5.4697 6.20267 (X¥X— 01307) = MOSS IFEX FX) 


SXELVEP GALLIUM SELENTOFE 
OINDEX(X) =SART (3.936242 0911 SEKEKS (XP K HM, BHB21% oSBB21) $10 PISUEKENS 
11K *X-160C.)) 
ETNDEXCXI=SORT C3. 313243 - SELB KMFEKS (XFK He 3H201% 23201) 410 TOT TEXOKS 
1(%*xX-1600.9) 


CADMIUM GERMANIUM ARSENTOF 
OLNNEXCKISSORT (Hs $90 AIL EK EK SUK EK He 55 26% 5524) +1 o BBBF KEK S (XK #X—1296. 
1)) 
FINDEX CX) ESORT (60 $90 S21 FKEKS (KEK = eG BG7* GHGT) 41 GDIEXF KS (X*K—1296.4 
1)) 


ZING GERMANIUM PHOSPHIDE 
OINDEX (X ESORT (9.74974 6 B89 (X* Ko 1402) OO 27REYFX) 
ETNDE XOX} =SORT (9.987 He 7U5LS (KKH 01489) = OO 271 KX) 


CINNARAR 
OTNDEX (CX IESART (CE. 944546 36587 (X* K-01 357) O01 IFEX FX) 
FINDEX(XISSART( 86 3922+ eS3ASCX* X01 38H) eo 0027 4X*X) 


PROUSTITE 
ATNDEX(XIESART(7 24822 + A3S/(X*X = 116) = 20016*X*X) 
FINDEX(XI=SQOT( Ae3434 + 2 33527(X#K = 21117) - -0007*X*X) 


PYRARGYRITE : 
OTNNEX (XI SSART (40 60 SASEXEKS(KX# KM o1G) He LAS SEXEKS (X*K—225.) ) 
ETNNEX(X)=SORT (10 $5 0 AUSEXFXS (MEK M16) He N20 2EKEKS (KF KH225,)) 


LIVHTUM NIOBATE 
OITNNEXCXI=HSORT(G.9264.117 SIX*#K 60473) — 0275 F XX) 
FINDEX(XISSORT (4657794. 9964S (¥*K=s0GKH) — oO 221 4K *X) 


THALLIUM ARSENIC SFLENTOE 
OTNOEXCXIASAAT (4A 0 H17 oA 2SEKEKS (KEK LUGS UGS) Hol FXFK/(K*KHI0—)) 
FUNDEXCXPSSART C10 t 8693 FXFKXS (KEK. GK5* 445) +o OSFXPK/S(KX*K—G000)) 


LITHIUM TODATE 
OTNDEX(XIESOPT (ZH 6 He O2ISS (KEK He DISA INH 77*XFX) 
FINDEXCXIESART( 2.961 te TL G27 (KEK 00915) — oe DO42TFKYK) 


CADMIUM SELENIDE 
ANTNDEX(XIESNOT(G, 1321410858 7EXENS IX EK. 2187) 3. OKbLEXEXS (KX FXHS3IAL.W 
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Cc 


SUBROUTINE DSQUAQ(DSQ, THETAR.I) 
COMMON/COEF /N(5) 


C THE FOLLOWING CARNS ARE FOR POINT GROUP &(BARD2M. THE COEFFICIENTS 
C MUST BE LISTED ON THE INPUT CARD IN THE ORDER 1(14)- 036). 


Cc 


Cc 


C THE FOLLOWING CAPDS ARE FOR POINT GROU? 32 


45 
20 


SINTHESIN(THETAR) 
SIN2Z2SIN(2.*THETAR) 
IF(I.NE.O) GO TO 10 
0$Q20(2)*0(2) ®*SINTH*SINTH 
RETURN 

IF(I<2) 15¢20925 
O5Q2.25*SIN2*SIN2Z* (D0 (1) 40(2)) #¥2 
RETURN 

0S02D (1) *0€1) *SIN2*SIN2 
RETURN 

95Q20(13*001) *SINTHESINTH 
RETURN 

ENO 


SUBROUTINE DSQUAR(DSO,THETAR,T) 
COMMON/CONEF /P (5) 


THE COEFFICIENTS MUST 


C PE LISTEN ON THE INPUT CARO IN THE ORNER (11)— 0(14). 


C 


SINTH=SIN(THETAR) 
COSSO=1.-SITNTH*SINTH 

TF(I.NE.0) GN TO 16 
DSQ=D0(1)*N(1)9¥*COSSO 

RETURN 

IF (T=#2) 15620425 
SIN2=SIN(2.*THETAR) 

ASN= (NCL) FC OSSH+ .5*N(2) FSINZ) ¥¥2 
RETURN 

SIN2Z=SIN(2.*THETAR) 
nsag=(0(1)*C0SS+ NC2V*SINZIEKS) 
RETUPN 

END 
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# 
SUBROUTINE OSQUAR (0S. THETAR,T) 
COMMON/COEF /0(5) 
Cc 
C THE FOLLOWING CARDS ARE FOR POINT GROUP 3M, THE COEFFICIENTS MUST 
C BE LISTED ON THE INPUT CARD IN THE ORNER (2104 N¢31).0015) 
C 
SINSQz (SIN(THETA®)) #¥2 
COS SQ2#1.-SINSQ 
SIN2=2SIN(2.* THETA) 
IFCI.NE.0) GO TO 10 
0SQ20(1)*0(1) ¥COSSO40 (204062) FSINSQ4N (1) #002) FS IN2 
RETURN 
10 IF(LT-2) 15415420 
15 O0SQ2COSSO*COSS0*0(1)*D¢1) 
RETURN 
20 DSQ20(3)4063) FSTNSO4N(1)90(1) FCOSSI4D (1) F003) *SINZ 
RETURN 
END 


SUBROUTINE ODSQUAR(NSO, THETAR,T) 
COMMON/COEF /N (5) 
Cc 
GC THE FOLLOWING CAPNS ARE FOR POINT GROUP 6 . THE COEFFICIENTS MUST 
C RE LISTEO ON THE INPUT CARO IN THE ORDER 0115), 0131), Di). 
C 


SINTHSSTI'C THE TAR) 
SINZ=SIN( 26 *THET AR) 
IF(INE*0) GO TO 10 
0sd=0(2)9*002) *SINTHFSINTH 
RETURN 

210 IFCT22) 15420425 

15 OSQ=.2570(3)*0(3)*SIN2*SIN2 
RETURN 

20 nsQ= NCZV*FN(SPFSINZ*SIN? 
RETURN 

25 0$0=0(1)*09¢€1) *STNTH*SINTH 
RETURN 
END 
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SURROUTTINE DSQUAR(OSQ.THETARsT) 
COMMON/COEF /N(5) 


c 
C THE FOLLOWING CAPNS ARE FOR POINT GROUP 6MM. THE COEFFICIENTS MUST 


C RE LISTED ON THE INPUT CARN IN THE ORNER 1(25)¢ 0NC31). 
Cc 
SINTHESIN(THETAP) 
IF(T.NE.0) GO TO 10 
NSQ#20(2)*0¢2) *SINTH*¥SINTH 
RETURN 
10 TF(I-2) 15.15.25 
15 OSQs6. 
RETURN 
25 OSQSD(1)*001) *SINTHFSINTH 
RETURN 
END 
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